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Abstract. In the present paper, we solve three boundary value problems related to 
the temperature field in oil strata - the fractional extensions of the incomplete lumped 
formulation and lumped formulation in the linear case and the fractional generalization 
of the incomplete lumped formulation in the radial case. By using the Caputo differin- 
tegral operator and the Laplace transform, the solutions are obtained in integral forms 
where the integrand is expressed in terms of the convolution of some auxiliary func- 
tions of Wright function type. A generalization of the Laplace transform convolution 
theorem, known as Efros' theorem is widely used. 
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1. Introduction 

The heat equation for a porous medium is derived in MUIOII under several assumptions 
on the modeL Besides the exact formulation of the temperature field in oil strata, the 
following three approximate formulations are also treated: 



• The lumped formulation. In this case, stratum's thermal conductivity is infinitely 
large in the vertical direction and is finite in the horizontal direction. The cap and 
base rock is considered to be thermally isotropic. 

• The incomplete lumped formulation. The horizontal heat transfer in both the cap and 
base rocks surrounding the strata is neglected. 

• The Lauwerier formulation. The horizontal heat transfer in the strata is also 
neglected. 

In the linear case of the incomplete lumped formulation the temperature field u — 
u{x,z,t) satisfies the differential equation given as follows (H]: 

— =a^^r-^, 0<x,z,t <o°,a>0, (1.1) 
dt dz 

subject to the boundary condition 

du d^u du du 
z = Q; ^-27— + a^, Q<x,t<oo^ (L2) 
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and the additional conditions 

(a) z = X = 0: M = 1 , 

(b) f = 0: M = 0, 

(c) M-^Oasjc^+z^^°o, (1.3) 

where 7 and a are positive constants. 

In the case of the lumped linear formulation, the temperature field problem is described 
by (1.1), (1.2) and (1.3) where (1.3)(a) is replaced by the following condition: 

z=x = 0; ^-lyu = -ly. (1.4) 
ax 

By using the Laplace and the Hankel transforms, similar boundary value problems with 
more general conditions were studied by Ben Nakhi and Kalla Some fractional gen- 
eralizations of the above discussed problem have been given by Boyadjiev and Scherer 
I3I4I . In the present paper we aim at further extending the works of Ben Nakhi, Kalla, 
Boyadjiev and Scherer. 

In this paper we shall consider the following fractional differential equation in the linear 
case of the incomplete lumped formulation u — u{x,z,t;l5), 

2p , d^u . du - 1 

M = fl^^-^ - A^: < jc,z < 00, < p < - (1.5) 
<3Z^ dz 2 

subject to the condition 

TR d^u du du 
z = 0:Dfu^^-2y— + a^, 0<x,f<oo, (1.6) 
ox^ dx dz 

and the additional conditions 

(a) z = x = 0: u— Ah{t), 

(b) if + — > 00, then u—^0, 

(c) f = 0:M = 0, (1.7) 

where A , A, 7, a are positive constants. 

We also discuss the linear case of the complete lumped formulation when the problem 
remains the same as given by (1.5) subject to the boundary conditions (1.6) and the con- 
ditions (1.7)(b) and (c) with (1.7)(a) replaced by 

z=x = 0:^-2Yu^-2yA. (1.8) 
dx 

We shall also study the following fractional incomplete lumped formulation in the radial 
case. Equation (1.5) in radial coordinates is considered as subject to the condition 

26 d^u 1 —2v du du 

z = 0:Dfu = ^ + — + a^, 0<r,f<oo, (1.9) 

dr^ r dr dz 
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and the additional conditions 

(a) z^r = 0:u=Ah{t), 

(b) if r +2" ^ then M — > 0, 

(c) t^O:u^O (1.10) 

Here denotes the Caputo fractional derivative f9\ defined as follows: 

{1 f' f"Hx) 
1^, ^ / \ I , d^: m-l < p <m 

K H> yt X) 
At'" 

The following rule of the Laplace transform will be used: 

;n-l 

L[D:/(0]=/l(/(0]-E/(0)/-^-\ m-l<^<m, (1.12) 

k=0 

where the Laplace transform of a function /(f) is defined as 

f{p)=L[f{t)]= f ^-P'f{t)At, Re(/7)>0. (1.13) 



2. Auxiliary results 

Efros' theorem [1|. If G{p) and q{p) are analytic functions and 

F{p)=L[f{t)],c-''^P^'G{p)=L[g{t,z)], (2.1) 

then 



F{q{p))G{p)^L 



/(T)g(f,T)dT 



where L is the Laplace transform operator 



(2.2) 



Efros' theorem is a generalization of the Laplace transform convolution theorem when 
q{p)=p. 

The following results will help in solving our problem 03181 : 

If < f,T,z < °o, 



and 



/(f,T;j3)=A^(^; jS), < j3 < i (2.4) 
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then 

(i) e-^p'' =L{g'it,T,zm, (2.5) 

(ii) ie-^'^"'=L{/(r,T;/3)}. (2.6) 

The above resuhs follow directly from p. 22, Lemma 2 of f3l| on replacing b = a/a. 
Here, functions M and are auxiliary functions defined as 

and 

= ^J ^e---'-'''d(T, 0<P<\, (2.8) 

where H denotes the Hankel path of integration that begins at (7 = — ibi{bi > 0), 
encircles the branch point that lies along the negative real axis and ends up at a = —0° + 

ib2{b2>0). 

The above functions are special cases of the following Wright function (or Bessel- 
Maitland function) lISTTTIl 

W{z;X,^i) = f .^.f ^. =^.1 ^^^dg, A>-1,M>0 (2.9) 

which is an entire function of z- 

From eqs (2.7), (2.8) and (2.9) the following relations can easily be deduced: 

M(z;j3) =W(~z;-j3,l-i3) (2.10) 

and 

A?(z;j3) =W(-z;-2/3,l). (2.11) 
For )3 = 1 /2, 

n(^z;^^ =Hit~z). (2.13) 

Here H denotes the Heaviside step function. 
For j8 = 1 , we have the known result IS) 

M(z,l)-5(z-l), 

where 5 is the Dirac delta function. 
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Lemma. IfO<t,T,z<°°,then 

(a) L-He-'""')^^M[-^;2I5], < 2^ < 1 (2.14) 

and 



•M(^-2^;2j3jdM, /orO<2j3<l. (2.15) 

Proof. Part (a) can easily be obtained on replacing )3 by 2)3 (0 < 2j3 < 1) and z — Q',a = 
a = 1 in (2.3) and (2.5). 

To prove part (b) we apply Efros' theorem. We first express 



_(z+«^),/„2/i^ll 

^- ^F{q{p))G{p)^ 



where G{p) = \,q{p) = p^^ + (AV4fl^) so that F{q{p)) = g-^^ 
Using a known result |[6l, we get 

F{p) = e « v'' = L — -=e. 4a2, 



Hence G(p)e-^i'?(P) = ^-T,(p^H().'-IAa^)) ^ -^he^^ 

.-'(G(ri.-*>)=e-^|iL«(iL:2?) 

By Efros' theorem we can easily get the result (2.15). 
In the case of 2)3 = 1, on using the formula lH) 

M(z,l)-5(z-l), 

where 5 is the Dirac delta function, we get, after some manipulations and applying the 
properties of the Dirac delta function, the following inverse Laplace transform 



L-|e-^V-^ Uli±^.-3/v^e-^. (2.16) 



3. Fractional extension for incomplete lumped formulation in the linear case 

Theorem 1. The solution of the linear case of the fractional incomplete lumped formu- 
lation problem (1.5), (1.6) and (1.7) is given by 



uix,z,t)^'j^J^ 2^__eS?(^+«^)gi(r,T,z;)3)dT, (3.1) 
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where g\{t,T,z;P) is defined as the Laplace transform convolution of three flmctions 
given by eq. (3.17). 

Proof. Taking the Laplace transform of eq. (1.5), boundary condition (1.6) and the con- 
ditions (1.7) with respect to variable t, we obtain by using the result (1.12), 

/'^M = a2 -A — :0<x, z<°o, (3.2) 
oz oz 

TO d^u ^ du du ^ 
z = 0: p^Pu = ^-2Y^ + a^,0<x<^, (3.3) 
dx-^ dx dz 

(a) z = x = 0: u^Ah{p), (3.4) 

(b) M^Oifx^+z^^oo, (3.5) 

where u = u{x,z,p;P) = L{u{x,z,t;P);t p} and h{p) denotes the Laplace transform 
of h{t) which occurs in the boundary condition (1.7)(a). The solution of (3.2) remains 
bounded as z — > oo and is given by 

u{x,z,p;l5)=B{x,p)c-i^, (3.6) 

where 



K = -^ + \ — + -^. (3.7) 
a^ V 

Using (3.6) in eq. (3.3) to determine function B{x,p) we get the ordinary differential 
equation 

under the conditions obtained from (3.4) and (3.5) as 

B{0,p)=Ah{p) and lmiB{x,p)=0. (3.9) 

The solution of (3.8) with (3.9) is given by 

B{x,p) =Ah{p)A'-^^'^^^^). (3.10) 

Thus the solution given by eq. (3.6) is expressed as 

u{x,z,p;l5) = Aei^^h{p)c-i'' - c-'^f+p'^+^K_ ^^j) 

We shall use Efros' theorem to find the inverse Laplace transform of (3.11) to obtain the 
solution of our proposed problem. 
Express 

u{x,z,p;li) = F{q{p;l5))G{p), 
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where 



Gip)=Ae^''h{p)e-i'^, 
qip)^f + p^P + ^K, 



= i{/(0}, 



(3.12) 



-^(''^/^)G(p) = Aer-r^^e-''^^(/,)e-^ 
= L(g(f, T,z;j3)). 



(3.13) 



Then 



u{x,z,t;l5)=L / f{T)g{t,T,z;l5)dt 

■'o 



(3.14) 



where functions / and g are defined in (3.12) and (3.13) respectively. 
Now to find g(f , T,z; j3), we write 



g(f,T,z;j3)=Ae^--'^^e^<''+"^'gi(f,T,z;j3), 



(3.15) 



where 



gi{t,T,z;li)^L- 



-rp- 



,2/i 



h{p)e 



(3.16) 



Using the Laplace transform convolution theorem. 



gi{t,T,z;l5)=L-\e-'p'^)*L-\h{p))*L-'ic''~'""''y'^' ' ^- 



(3.17) 



Here in (3.17) the functions can be simplified by using results (2.14), (1.7)(a) and (2.15). 
Hence the required solution is given as 



u{x,z,t;j5) = 



e^'^'+"''gi(f,T,z;j3)dT. 



(3.18) 



Special cases 

(i) When 2j3 = 1, Theorem 1 reduces to the problem studied by Ben Nakhi and Kalla 121 
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and the solution can be obtained by representing the function gi in (3. 17) as 



jz+ar) 



^ {5(f-T) */!(?)}* 

= {h{t-z)H{t-z)} 

_ H{t-x){z + ax) 
2a\f% 

On substitution of (3.19) in (3.18), we get 



-i — ^'^e 4a^' e 

2flY'7r 



_3/2 -(£±££11 -4, 

' e 4fl2f g 4a2 



u ^' e 4fl'2» e 4?"/!(f— T— M)dM. 

(3.19) 



Aan 



r3/2 



-e2a- 



r(z+at) 



dr. 



(3.20) 

(ii) On taking A = 0,A = /z(f ) = 1 in (1.5) and (1.7)(a) Theorem 1 reduces to the problem 
considered by Boyadjiev and Scherer f3 | and the solution is obtained by representing the 
function as a convolution of two functions in Theorem 1 given as 



gl{t,T,Z\^)=L 



/ rn-P 



From (2.5) and (2.6), we get 

gi(f,T;z;j3)=g'(f,T,z;j3)*/(f,T,z;j3). 
If we further take /3 = 1/2, it yields the result given in Antimirov et alfV\. 



(3.21) 



(3.22) 



4. Fractional lumped formulation in the linear case 

Theorem 2. The solution of the fractional lumped linear formulation (1.5), (1.6), 
(1.7)(b), (c) and (1.8) is given by 

•e^*^'+"'*g2(f,T,z;j3)dT, (4.1) 

where g2{t, is defined as the Laplace transform convolution of two functions given 

by (2.6) and (2.14). 
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Proof. Taking the Laplace transform of the equations of the fractional lumped formula- 
tion problem, it reduces to the following boundary value problem: 

p'^Pu = a-^^r-^-X — , 0<x,z<°°, (4.2) 
oz oz 

„ 9fi d^u ^ du du ^ 

z = 0;p^Pu = ^-ly^ + a^; 0<jc<oo, (4.3) 
ax'^ ax az 

du ly 

(a) z = x = 0: ^-2yu = -A, 

ax p 

(b) ifx^ + z^^oothenS^O. (4.4) 
Proceeding as in Theorem 1, the solution to (4.2) is 

u{x,z,p;P)^B{x,p)c-i^, (4.5) 

where 



K=-A + -\P'^ + ^- (4.6) 
2a^ a y 4a^ 

It is bounded as z — > «>. 

Now at z = 0, (4.5) satisfies eq. (4.3) i.e. 

(fB . dB 



— T - 27 

dx^ dx 



2\ 



B = 0. (4.7) 
As limx^ooB(x,p) = and B satisfies (3.8), solution to (4.7) is 

B{x,p) =B{p)A'-^^^^^),B{p)=B{0,p). (4.8) 
From (4.4)(a) and (4.8), 

B{p) = — . . (4.9) 

Hence from (4.8), 

Bix,p) = 2yA ^jr-VfT?^) 

and from (4.5), it is expressed as 

uix,z,p;P) = - Ar-V^^^),-iK_ (4.11) 

p[Y+^y^+P^^ + ^K) 

Now using Efros' theorem to express u{x,z,p\fi) — F{q{p;j5 j)G{p), we take 

q{p;^) = y'+p'^ + ^K, 
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and 



F{q{p;li)) = 



Using formula ||6l 

F{p) = 



Now 



y+^/p 



G{p;j3) = 2yAe^'-'^'e^ 



±^,1 pip + 22, 

" V " 



(4.12) 



where 



ife(f,T,z;/3)] 



.2 if74 



G(/,;j3) = 27Ae''---r e5:?^-'^"^'L[g2(f ,T,z;i3)], 



-1 



r-1 



-T/r 



,2)3 - 



-1 



(4.13) 
(4.14) 



(4.15) 



Functions in the right-hand side of (4.15) are represented by (2.6) and (2.15) respectively. 
Thus the solution (4.1) is obtained by using results (4.12), (4.13), (4.14) and (4.15) in the 
following: 



M(x,z,f;j3) = / /(T)^2(f,T,z;j3)di 

•J n 



Special cases 

(i) When 2j3 = 1, Theorem 2 reduces to the form studied by Ben Nakhi and Kalla ijS). To 
derive it we take 2j3 = 1 in (4.15) and use the result (2.16) to get 

= i/(r-T)Jii±^rVV^e-Sl 
^ ' \ 2aVn J 

zj^ ^ + -3/2 - '-^^ 
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Now using the integral given in HI we get 

g2 ^ // f-T e^?^ 'erfc^^ ^ 

+ e 'erfc ^ — ; - — yf- T . 

Substitution of giit^ j3) in (4.1) gives us the desired result. 

(ii) For A = 0,A = h{t) = 1 and setting b^a/a the problem (1.5), (1.6), (1.7)(b), (c), 
(1.8) and the corresponding solution given by (4. 1 ) reduce to the fractional formulation in 
the lumped linear case studied in (O . 

Further, in case j3 = 1/2, our problem reduces to the form considered in Antimirov 
et al yj. 



5. Fractional incomplete lumped formulation: The radial case 

Theorems. Here we consider temperature field problem for fractional incomplete 
lumped formulation in the radial case represented by (1.5), (1.9) and (1.10) and the 
solution is given by 

u{r,z,t-l5)^-^ [ -l-fn''e-5^e^(^+"^V(f,T,z;i3)dT, (5.1) 



T{v)J, tv+1 V2 

where gi(f,T,z;j3) is the Laplace transform convolution of the three functions given by 
eq. (3.17). 

Proof. Proceeding as in the previous theorem, if we take the Laplace transform of prob- 
lem (1.5) and the set of conditions (1.9) and (1.10), we get the following: 

p^Pu^a^^~X — , 0<r, z<oo, (5.2) 
oz^ oz 

90 (9^M \ — 2v du du 

z = Q:p^Pu = ^^ — + a— , 0<r<oo, (5.3) 

ar^ r or az 

(a) r = z = Q:u=Ah{p), 

(b) if r2+z^-^°°thenM-^0. (5.4) 
The solution of (5.2) which remains bounded as z — > 00 is 



«(r,z,p;i3) = B{r,p)^-^^,K - + -\ P^^ + (5-5) 

a^ ay Aa^ 

To determine the function B{r,p), substitution of (5.5) into (5.3) gives the following ordi- 
nary differential equation: 



a 



d^B l-2vdB 
dr^ r dr V 2 



p^l^ + ^K]B = Q (5.6) 
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and the corresponding initial conditions are obtained by using (5.5) in (5.4) as follows: 

B{0,p)=Ah{p) and lim B(r,p) = 0. (5.7) 

Now the bounded solution of eq. (5.6) with the condition (5.7), using the formula 
(p. 986, 8.494 (9) of |7|) is given as 



(5.8) 



where ji = ^ p^P + {a/2)K and Ky{jj.r) is the modified Bessel function of the second 
kind, of order v, and Bi (p) is a constant which can be determined with the help of (5.7) 
and the standard result. 



^vW~^r(v)0) asz^O 



(5.9) 



IS given as 



Blip) 



r(v) 



Ahip). 



Now substitution of (5.10) and (5.8) into (5.5) gives the solution: 
M(r,z,p;j3) = :^Ah{p) (^) Vv(Mr)e-f . 



(5.10) 



(5.11) 



To get the solution of our proposed problem, we shall find the inverse Laplace transform 
of (5.11) with the help of Efros' theorem. So, let us first express 



where 



and 



a{r,z,p;l5) = G{p-l5)F{q{p-l3)), 



G{p;fi) 
F{q{p-m 



r(v) 



Ah{p)e 



-iK 



fir 



(5.12) 



(5.13) 



^^y/qip;p). 
Using the well-known result ||6l, we get 



Fip) 



2 

mm 



Ky{r^7)=L 



1 /r\2v 1 _,2 
2V2/ 



(5.14) 



Fractional extensions of some boundary value problems in oil strata 279 



Now 

2A 



= L 



2A A, 



r(v) 

= L[g{t,T,z;p]. (5.15) 

Then by Efros' theorem, the solution (5.1) is obtained by using (5.14) and (5.15) into the 
following 

u{r,z,t;l5)= [ /(T;/3)^i(;,T,z;/3)dT. (5.16) 

Here ^i(/,T,z;/3) is the Laplace transform convolution of the three functions given by 
eq. (3.17). 

Special cases 
(i) 

COROLLARY 1. 

If we take 2/3 = 1 in Theorem 3, the problem reduces to the following form: 

^=a%^-A^; 0<r,z,?<oo, (5.17) 
ot az az 

du d^u l—2vdu du 

z = 0:— = ^H — + a^; 0<r,?<oo, (5.18) 

at dr'^ r dr dz 

(a) r = z = 0: M = A/i(?), 

(b) if + ^ oo, then m ^ 0, 

(c) ; = 0:m = 0, (5.19) 
and the solution is given by 

2 

M(r,z,r;l/2) = ^^ (^^j i^e^(^+«^)gi(^T,z;l/2)dT, (5.20) 

where 



__(j+aT) / a2 

-1 o-Tp/,c^^o 4? 



gi(^?,T;z;i)=L-i|^e-^/';j(;,)e 

= L-i(e-^''ft(p))*L^' ( e " 
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Using 

L-\e-^Ph{p))^h{t^'z)H{t^x) (5.21) 
and (2.16) the solution is expressed as 

, 1\ A / 1 (r^y _^ Mz+ar) 



/ ^ ■j=^u ' h[t-T-u)e 4^e ^ATdit dr. 



(5.22) 

(ii) On taking A=0,A=/i(f) = l and setting = a/a in Theorem 3, we obtain the result 
given by Boyadjiev and Scherer |f3l. 

Further for 2)3 = 1, we get the result given in Antimirov ef [T]. 
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